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The production of strange pentaquark states (e.g. 
Theta baryons and S states) in hadronic interactions 
within a Gribov-Regge approach is explored. In this ap- 
proach the 0"'"(154O) and the S are produced by disinte- 
gration of remnants formed by the exchange of pomerons 
between the two protons. We predict the rapidity and 
transverse momentum distributions as well as the 47r mul- 
tiplicity of the e+, S— , S", 2° and S+ for = 17 GeV 
(SPS) and 200 GeV (RHIC). For both energies more than 
10^'^ 8+ and more than 10^^ S per pp event should be 
observed by the present experiments. 

Very recently in photon-nucleus [1,2] and Kaon-nucleus 
experiments [3] a new baryon, consisting of five quarks, 
uudds, has been identified in the K^n or K^p invariant 
mass spectrum. It has been been named particle and 
has spin 1/2, isospin and strangeness -1-1. Its mass 
is about 1.54 GeV and its width is less than 25 MeV. 
Such a state has been predicted by Diakonov [4] in the 
framework of a chiral soliton model. 

This finding has renewed the experimental and theo- 
retical interest for novel baryon states. Major progress 
has been reported on a possible extension of the original 
model [5,6], within the Skyrme model [7-9], and within 
the constituent quark model [10,11]. Also (lattice) QCD 
studies of the 0+ (see e.g. [12-14]) have been performed 
and first explorations of the 8+ multiplicity at SPS and 
RHIC energies are available [15-17]. Thus, the existence 
of this novel state has many perspectives for pp as well 
as for nucleus- nucleus collisions. 

In this letter, we present predictions for the 8+ in pp 
collisions from a newly developed approach for hadronic 
interactions. It has recently been shown [18] that the 
standard string fragmentation models which described 
spectra and multiplicities of many hadrons rather well 
need to be revised: due to their diquark-quark topol- 
ogy these models produce more Q than f2 in medium 
and low energetic pp interactions, in contradistinction to 
experiments. Therefore a key issue is presently to gain 
information on the details of hadron production. Espe- 
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cially more exotic states like the A(1405) which may be 
a udsuu state or the 8"*" carry important information to 
tackle this question. If their multiplicity can be related to 
that of other particles one can hope to get experimentally 
a handle on the hadronization process. 

In heavy ion collisions the multiplicity of the most 
abundant particles can be well described in a statisti- 
cal model assuming a temperature close to that where 
the chiral/confinement phase transition is expected and 
a moderate chemical potential. Unstable particles can 
test how the expanding system interacts afterwards be- 
cause if the decay products have still an interaction in the 
invariant mass spectra the resonance cannot be identified 
anymore. Especially long living states are very useful in 
this respect. 

Until this year, the search for multi-quark bags focused 
mainly on the _ff-particle [19], because it is closely related 
to the study of S and AA hypernuclei (see e.g. [20,21]). 
The H is a. six quark state (uuddss) coupled to an SU(3) 
singlet in color and fiavour. Unfortunately no stringent 
observation of the iJ-particle exists. Even today, decades 
after the first prediction of the 77-di-baryon by Jaffe [19] 
the question of its existence is still open. 

In contrast to the H particle, the situation for the 8+ 
baryon is very promising. Thus, in this letter we explore 
the formation of the 8"'"-baryon within a new approach 
called parton-based Gribov-Regge theory. It is realized 
in the Monte Carlo program NeXuS 3.97 [22,23]. In this 
model high energy hadronic and nuclear collisions are 
treated within a self-consistent quantum mechanical mul- 
tiple scattering formalism. Elementary interactions, hap- 
pening in parallel, correspond to underlying microscopic 
(predominantly soft) parton cascades and are described 
effectively as phenomenological soft Pomeron exchanges. 
A Pomeron can be seen as layers of a (soft) parton ladder, 
which is attached to projectile and target nucleons via leg 
partons. At high energies one accounts also for the con- 
tribution of perturbative (highpi) partons described by a 
so-called "semihard Pomeron" - a piece of the QCD par- 
ton ladder sandwiched between two soft Pomerons which 
are connected to the projectile and to the target in the 
usual way. The spectator partons of both projectile and 
target nucleons, left after Pomeron emissions, form nu- 
cleon remnants. The legs of the Pomerons form color 
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singlets, such as q-q, q-qq or q-qq. The probabihty of 
q-qq and q- qq is controlled by the parameter Pqq and 
is fixed by the experimental yields on (multi-) strange 
baryons [23]. 
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FIG. 1. The typical collision configuration has two rem- 
nants and one cut Pomeron represented by two q — q strings. 
In the case of a uudds flavour content a 0"*" baryon can form. 



Particles are then produced from cutting the Pomerons 
and the decay of the remnants. As an intuitive way to 
understand particle production, each cut Pomeron is re- 
garded as two strings, i.e. two layers of a parton lad- 
der. Each string has two ends which are quark(s) or 
antiquark(s) from the two Pomeron legs respectively. To 
compensate the flavour, whenever a quark or an anti- 
quark is taken as a string end, a corresponding anti- 
particle is put in the remnant nearby. 

Since an arbitrary number of Pomerons may be in- 
volved, it is natural to take quarks and antiquarks from 
the sea as the string ends. In order to describe the ex- 
perimental yields on (multi-)strange baryons [23], all the 
valence quarks stay in the remnants, whereas the string 
ends are represented by sea quarks. Thus, Pomerons 
are vacuum excitations and produce particles and anti- 
particles equally. Note that in addition to these sin- 
glet type processes, valence quark hard interactions are 
treated differently in the present model. To give a proper 
description of deep inelastic scattering data, a certain 
fraction of the Pomerons is connected to the valance 
quarks of the hadron, not leading to a quark feeding of 
the remnant. This kind of hard processes will not be 
discussed here, but is included in the simulation. Only 
the remnants change the balance of particles and anti- 
particles, due to the valence quarks inside resulting in 
the possibility to solve the anti-omega puzzle [18] at the 
SPS. 

This prescription is able to accumulate multiple quarks 
and anti-quarks in the remnants depending on the num- 
ber of exchanged Pomerons. In the most simple case of a 
single Pomeron exchange, the remnant may gain an ad- 
ditional anti-quark and a quark and is transformed into 
a pentaquark bag as discussed in the following. 

The typical collision configuration has two remnants 
and one cut Pomeron represented by two q — q strings, 
see Fig. 1. Four quarks and one anti-quark are now in 
the remnant: the three valence quarks u, u, d plus a q^q^ 
sea quark pair. Each of the q^ may have the flavour u , d, 



or s, with relative weights 1 : 1 : /s. The parameter /s, is 
fixed by strange hadron data in proton-proton scatterings 
at 160 GeV to be fs = 0.26. Thus, there is a small but 
nonzero probability to have a uudds fiavour in a remnant, 
such that a 0+ pentaquark may be formed. 
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FIG. 2. Rapidity distributions of 0"'"'s in pp interactions 
at Ejg^j^ = 160 GeV. 
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FIG. 3. Transverse momentum distributions of 0^'s in pp 



interactions at E- 
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160 GeV. 
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remnants have mass distribution P(m^) cx 
, TO^ e (Tij^nin' ^^'^rc s IS the squared CMS 

energy. With, TOmin being the minimal hadron mass com- 
patible with the remnant's quark content, and is the 
light-cone momentum fraction of the remnant which is 
determined in the collision configuration. In the present 
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study, the parameter a is 2.25 for non difFractive inter- 
actions and 1 for diffractivc events [24]. This remnant 
disintegrates into hadrons according to (microcanonical) 
phase space [25]. This approach describes quite well 
multiplicity^ , transverse momenta and rapidity of all the 
observed hadrons [24] in pp collisions at energies between 
40 and 160 GeV. 

For the present study we have embedded the 0+ in the 
microcanonical approach to describe the disintegration 
of the remnant. It is therefore treated the same way as 
all the other hadrons. A similar approach was recently 
introduced to study H° production in pp interactions 
[26]. 

As discussed above, a pp collision in the present model 
involves two sources of particles production: 

1. Fragmentation of strings (from exchanged 

Pomerons). 

2. Decay of the projectile/target remnant. The decay 
of the remnant can happen in two distinct ways: 

2a.) by string fragmentation or 
2b.) by micro-canonical phase space. 

Presently, Pentaquarks can only be produced from 

the micro-canonical remnant decay. A fragmentation of 
strings into Pentaquarks is not included. 

This is in contrast to other hadrons, e.g. Lambdas that 
can be produced from all mentioned sources. 

There is no difference in the general input parame- 
ters for calculations for As and Pentaquarks. The main 
parameter is the probability for having a (anti-)strange 
quark at the Pomeron leg. This parameter controls 
strangess (e.g. Lambda and Kaon multiplicities) and 
Pentaquark production directly. 

The additional inputs for the study of Pentaquark 
states are their masses, quark contents and spin, which 
have ultimately to be determined by the experimental 
results. Presently, we use the values from Ref . [4 - 6] . 

Contrary to grand canonical thermal model estimates 
[15,17], which implicitly assume chemical and thermal 
equilibrium, the present approach works differently: In 
the present approach, initial and final state dynamics of 
the pp interaction are fully taken into account, thus al- 
lowing to predict transverse and longitudinal momentum 
distributions all produced hadrons, including the 0+. 
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FIG. 4. Rapidity distributions of 0"^'s in pp interactions 
at ^/s = 200 GeV. 
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FIG. 5. Transverse momentum distributions of Q'^'a in pp 
interactions at -^/s = 200 GeV. 



^A variation of the model parameters by 20% results in vari- 
ation of the model results by roughly 10% both for the Pen- 
taquarks and the well known hadrons. A larger variation of 
the parameters destroys the agreement with the NA49 pp data 
at 160 GeV. 



Particle 


Quark content 


Yield (SPS) 


Yield (RHIC) 


0+ 


(uudds) 


3.5 X IQ-" 


5.2 X IQ-" 




(ssuud) 


4.8 X IQ-^ 


11.9 X IQ-^ 




(ssudd) 


4.7 X IQ-^ 


9.2 X IQ-^ 




(ssudu) 


3.0 X lQ-5 


6.3 X IQ-s 




(ssddu) 


1.8 X IQ-^ 


4.5 X IQ-^ 



TABLE I. Predictions of the various Pentaquark abun- 
dances in 4-K for inelastic pp collisions at i?iab = 160 GeV 
(SPS) and = 200 GeV (RHIC). 
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Furthermore, a microcanonical model is employed for 
the remnant decay which counts exactly the phase space 
states for a given energy and volume of the system. Only 
in the limit of large volumes and many particles it agrees 
with the grand canonical approaches. 

Let us now study the multiplicities and momentum 
spectraof the calculated 8+ 's. Fig. 2 depicts the rapidity 
distribution of the predicted 6+'s at the top SPS energy. 
One observes a slight dip around central rapidities due to 
the fact that the particles originate from remnants. The 
total 6+ multiplicity per inelastic event is 0.0035. Table I 
gives also the predicted yields of the other species of Pen- 
taquarks. It is interesting to note that the S Pentaquark 
is suppressed by two orders of magnitude compared to 
the Theta particle. This is due to the double strange 
quark component of the S Pentaquark that is suppressed 
at Pomeron legs. However, it is interesting to note that 
the double-strange S Pentaquark has been well observed 
at this energy by the NA49 collaboration. Thus, with the 
presently accumulated statistics of the NA49 collabora- 
tion, the 0+ should be easily visible in the data. Fig. 3 
depicts the transverse momentum spectra of the O+'s. In 
Figures 4 and 5 we show the rapidity and transverse mo- 
mentum spectra at ^/s = 200 GeV. Especially at RHIC 
energies one clearly observes the pile-up of baryons 
in the forward and backward hemisphere. In the midra- 
pidity region the pentaquark yield vanishes. It should 
be noted that the vanishing yield at central rapidities 
might be filled up, if also the (yet unknown) fragmen- 
tation function into Pentaquarks in a string break up is 
included. 

In conclusion, we have studied the production of 
Pentaquark baryons in pp collisions from parton-based 
Gribov-Regge theory. Multiplicities, rapidity and trans- 
verse momentum spectra for the 0+ and S Pentaquarks 
are predicted for pp interaction at Ejg^^^ = 160 GeV and 
= 200 GeV. We predict a total yield of 3 - 5 x 10"^ 
0+ per inelastic pp event at SPS and RHIC. The S pen- 
taquark states are suppressed by two orders of magnitude 
compared to the 0+. Thus, our predictions are accessi- 
ble by the NA49 experiment at CERN and the STAR 
experiment at RHIC. 
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